Finite size effect in the specific heat of YBa 2 Cu 3 07_5 
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We analyze high precision, high resolution microcalorimetric specific heat data of detwinned 
YBa2Cu307_,5 single crystals. It is shown that a finite size effect must be taken into account when 
attempting to extract critical properties. Contrary to previous claims, there is no contradiction to a 
regular second order phase transition. Consistency with 3D-XY critical behavior is demonstrated. 
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Recently, detailed high precision, high resolution 
microcalorimetric specific heat data of detwinned 
YBa 2 Cu307_,5 single crystals were reported by Char- 
alambous et al. Q. Neglecting corrections to scaling 
and inhomogeneity induced finite size effects it was ar- 
gued that these measurements show definite incompat- 
ibilities with 3D-XY critical behavior and provide evi- 
dence that the specific heat critical exponents are asym- 
metric around T c . Furthermore it was concluded that 
the measurements are incompatible with a regular sec- 
ond order phase transition. This claims contradict the 
hitherto well established evidence, that in the experimen- 
tally accessible temperature regime the critical properties 
of cuprate superconductors fall into the 3D-XY univer- 
sality class . 

In this letter we reanalyze the data of jfj] by taking 
an inhomogeneity induced finite size effect into account. 
We identify a reduced temperature regime bounded by 

tf < t+ i5 tph and C ~ * _ ~ tph- As usual, in tne 
above expression the ±-signs refer to t > and t < 0, 
respectively, where t = (T — T c )/T c . Below tf a micro- 
scopic inhomogeneity mediated finite size effect is iden- 
tified. It makes it impossible to enter the asymptotic 
scaling regime, while above t h the temperature varia- 
tion of the background, usually attributed to phonons, 
becomes significant. To set the scale we note that in the 
A-transition of 4 He there are no strains and the only sig- 
nificant impurity is 3 He, which occurs naturally only to 
the extent of 1 part in 10 7 , and which at that level has a 
negligible effect O]. Indeed, the critical properties can 
be probed down to t — 1CP 9 . In solids and, in par- 
ticular, in cuprate superconductors such a homogeneity 
can hardly be achieved. Even the best solid materials 
exhibit phase transitions which are rounded over a range 
of < 10~ 4 . For this reason finite size effects must be 
taken into account when an attempt is made to extract 
critical properties from experimental data. Our analysis 
clearly uncovers a finite size effect in the specific heat 
data of detwinned YBa2Cu307_^ single crystals El. On 
this basis it is shown that in an intermediate tempera- 



ture regime the data are fully consistent with 3D-XY 
critical behavior and that the neglect of the finite size ef- 
fect leads to erroneous conclusions concerning the critical 
properties. 

The letter is organized as follows: we sketch the scal- 
ing form of the specific heat including the finite - size 
effect, as well as the corrections to the leading behavior. 
On this basis we analyze the experimental data, identify a 
finite size effect in terms of a scaling function and demon- 
strate that in the experimentally accessible temperature 
regime around the transition temperature T c , limited by 
the finite size effect, there are no incompatibilities with 
3D-XY universality. 

With the exception of the A-transition in 4 He it is in 
practice very difficult to access the asymptotic critical 
regime t — > fll|] . Due to inhomogeneities, a solid is ho- 
mogeneous over a finite length L only. In this case, the 
actual correlation length = £o = |£|~ iy cannot grow 

beyond L as t — > 0, and the transition appears rounded. 
Because of this finite size effect, the specific heat peak 
occurs at a temperature Tp shifted by an amount 
compared to the homogeneous system. The magnitude 
of the peak located at temperature Tp scales as L 01 '" 
|l2|| . For finite L the specific heat coefficient adopts the 
scaling form 

C/T = C(t P ,L = oo) + L a ' v TM, y = tL 1 /\ (1) 

where T\ (y) is a universal function of its argument and 

C(t P ,L = oo) = B + (0,oo)+A+\t P f a . (2) 

The study of finite size effects in liquid helium near the 
A point, as the bulk is confined more and more tightly 
in one or more dimensions, has been a very good ground 
for testing the finite size scaling theory ]l3|-[lo]|. From 
this work, the form of the scaling function T\(y) is rea- 
sonably well known for various confining dimensions and 
boundary conditions. 

In homogeneous systems and near the critical point, 
the general from of the specific heat along path t can be 
written as |l6| 
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f s is the free energy density, and a are the lead- 
ing specific heat critical amplitude and exponent, re- 
spectively. > denotes the leading correction ex- 
ponent, ejr and B^ are the amplitudes of the cor- 
rection terms. These dimensionless correction amplitudes 
are nonuniversal. They determine the size of the critical 
region and will in general become larger near a crossover 
p7[ . Nevertheless, the leading behavior as t — > is \t\~ a , 
but even in the A-transition of 4 He this corrections must 
be taken into account when attempting to extract reliable 
estimates for the critical exponents and amplitudes. 
The specific heat coefficient and C s are related by 
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where F is the free energy and V the volume. When the 
hyperscaling relation, Zv = 2 — a, holds, various criti- 
cal amplitude combinations adopt universal values 0,ff6| . 
Examples include a+/a~ and 
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is the critical amplitude of the correlation length 
along i — (x,y,z). In Tab. [I] we list recent estimates 
for critical exponents and amplitude combinations of the 
3D-XY universality class. 

To clarify the relevance of an inhomogeneity mediated 
finite size effect in cuprate superconductors, we calcu- 
lated the scaling function J-~i(y) (see Eq. dl])) from the 
reported data Q|. The result shown in Fig. |l| exhibits 
the characteristic form of the scaling function for a sys- 
tem in confined dimensions, i.e. corresponding to rod or 
cube shaped inhomogeneities [Ej). As a consequence, our 
analysis clearly reveals that the YBa2Cu307_,5 samples 
considered here exhibit a finite size effect on a length 
scale L ranging from 300 to 400 A (Tab. ||). For this 
reason, deviations from 3D-XY critical behavior around 
t p (Tab. 0) do not signal the failure of 3D-XY univer- 
sality, as previously claimed jjj, but reflect a finite-size 
effect at work. Next we sketch our data analysis and the 
derivation of the estimates listed in Tab. 0. 

For this purpose we explore the occurrence of 3D-XY 
critical behavior in the intermediate temperature regime 
(i.e., away from tp ). Since a is small, (A^ / a)\t\~ a rj 
A ± /a — A ± ln\t\, it is advisable to plot the measured 
heat coefficient data versus In |t| or log 10 It). Such a plot 
is shown in Fig. |^ for sample YBC03 111. The upper 
one corresponds to T < T c and the lower one to T > T c . 
The filled circles in Fig. |^ mark the intermediate scaling 
regime, while the open circles closer to T c belong to the 
finite size affected regime. T c of the fictitious homoge- 
nous system was adjusted to yield optimal consistency 



with nearly straight branches in an intermediate regime, 
subjected to the finite size constraint T c > T p . In this in- 
termediate temperature region (filled symbols), bounded 
by the region where the finite size effect dominates and 
the regime where the phonon contribution leads to signifi- 
cant corrections. Given this T c , the remaining parameters 
A ± and B ± , entering the leading 3D-XY form 



C*/T = A ± 10" Qlogl « 1 * 1 +B d 



(6) 



were the fixed by a linear fit to the corresponding 
branches for a = -0.013 and A+/A~ = 1.07 (Tab. fl). 
The resulting estimates are listed in Tab. [n]. 

From Fig. ^| it is seen that this procedure provides 
a reasonable description and with that meaningful pa- 
rameters needed to determine the scaling function T\ (y) , 
depicted in Fig. [j]. Indeed, the resulting critical am- 
plitudes A^ are consistent with the previous estimate 
@, A- w 1.5m J I (gK), for a YBa 2 Cu 3 07-5 single crys- 
tal with T c m 93.05K. Further away from T c the de- 
viations can be attributed to a temperature dependent 
background contribution, accounting for nearly 95% of 
the specific heat coefficient. To illustrate this point we 
depicted in Fig. ^| the temperature dependence of the 
specific heat coefficient of sample YBC03 in a wider 
temperature range around T c . Hence, due to the finite 
size effect and the temperature dependence of the back- 
ground the intermediate regime is bounded by the tem- 
perature region where the data depicted in Fig. |^ fall 
nearly on straight lines. Clearly, this interval is too small 
to estimate the corrections to scaling (Eq. (|^)) reliably. 
In this context it should also be kept in mind that any 
subtraction of the temperature dependent background in- 
terferes with the corrections to scaling and that the es- 
sential limitation of the intermediate regime arises from 
the finite size effect. 

Given the estimates of the parameters listed in Tab. || 
it is now straightforward to construct the scaling func- 
tion J-\(y) and to calculate quantities of interest from 
the experimental data. Of particular interest is the char- 
acteristic length of the inhomogeneities. Noting that at 
T„ 
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holds, we can express L with aid of the universal ampli- 
tude combination (|^) in terms of A + and tp: 
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Interestingly enough, the resulting estimates, ranging 
from 300 to 400 A (Tab. ||) clearly point to rather mi- 
croscopic inhomogeneities. 

To summarize, we have shown that the specific heat 
measurements considered here are remarkably consistent 
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witt 3D-XY critical behavior. An inhomogeneity in- 
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duced finite size effect makes it impossible, however, to 
enter the asymptotic scaling regime, where the leading 
scaling term clearly dominates. Nevertheless, the ex- 
posed consistency of the specific heat data with 3D- 
XY critical behavior in the intermediate temperature 
regime, combined with evidence emerging from pene- 
tration depth, other specific heat and magnetic torque 
measurements, etc. |p|-|lO|, clearly reveals that the pre- 
vious claim JlJ: 'evidence for asymmetric critical expo- 
nents, violating 3D-XY universality', based on the iden- 
tical measurements, must be attributed to an inappro- 
priate data analysis, neglecting finite size effects. Con- 
sequently, finite size effects must be taken into account 
when attempting to extract critical properties from the 
data |P1P^|. Finally, the resulting limitation renders the 
experimental search for charged fixed point behavior in 
cuprate superconductors into an academic issue. 

The authors are grateful to P. Gandit for providing 
the data files and to J. Hofer, A. Junod and C. Meingast 
for very useful comments and suggestions on the subject 
matter. 
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TABLE I. Estimates of critical exponent and universal am- 
plitude combinations for the 3D-XY universality class. 





YBC03 


UBC2 


UBC1 


A+ 


\mJ/{gK 2 )] 


-1.52 


-1.52 


-1.52 


A+ 


[mJ/{gK)) 


1.820 


1.824 


1.822 


A+ 


[10- 20 cm 3 ] 


8.4 


8.4 


8.4 


B+ 


[mJ/(gK 2 )\ 


3.554 


3.562 


3.559 


E>~ 


[mJ/(gK 2 )} 


3.528 


3.544 


3.54 




Tc [K] 


92.12 


92.33 


92.29 




T P [K] 


91.98 


92.24 


92.21 




t P [10 3 ] 


-1.52 


-0.94 


-0.87 




L [A] 


290 


397 


419 



TABLE II. Estimates for the critical amplitudes A + , 
the background terms B , the transition temperature T c 
of the ficticious homogeneous system, the temperature Tp, 
where the specific heat coefficient reaches its maximum 
value, and the characteristic length L of the inhomogeneities. 



(A+ = T c A+a, 
p = 6.37#/cm 3 .) 



A+ [cm 3 ] = p/k B W 4 A+ \mJ/{gK)], 
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FIG. 1. Scaling function TxitL 1 /") derived from the exper- 
imental data of sample YBC03 (o), UBC2 (□) and UBC1 (o) 
with the aid of Eqs. ([!]) and (|^) using the parameters listed 
in Tab. [fl|. 
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FIG. 2. Specific heat coefficient C/T [mJ/(gK 2 )\ ver- 
sus log 10 \t\ for YBa2Cu307_« (sample YBC03) with 
T c = 92.12^. The filled symbols denote the 3D-XY scal- 
ing regime, whereas the open circles belong to the finite size 
and background affected regimes, respectively. The solid lines 
are fits to Eq. (H) using the filled data points only. 
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